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Abstract 

Collective low lying levels of light and medium Xenon isotopes are deduced from the Generalized Bohr Hamiltonian (GBH). 
The microscopic seven functions entering into the GBH are built from a deformed mean field of the Woods-Saxon type. 
Theoretical spectra are found to be close to the ones of the experimental data taking into account that the calculations are 
completely microscopic, that is to say, without any fitting of parameters. 
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I. INTRODUCTION: 



The so-called General Bohr Hamiltonian (GBH) is applied to the light-medium even-even Xenon isotopes region 
112 < A < 126. This region lies between very deficient neutron nuclei with half-lives of few seconds and stable nuclei 
(A = 124, 126). 

The lightest isotopes (A = 112 — 116) are not very far away from the drip line (exotic nuclei) and hence experimental 
data such as band sequences and probability transitions are difficult to obtain. For these nuclei, we can find some 
data in Ref. Data of other isotopes (A = 120 - 126) are given in Ref.[l]-[l||. 

Some simple properties can early be deduced from the concept of the equilibrium deformation. In effect, these isotopes 
have a proton number (Z = 54) which is somewhat close to the magic number of a closed shell (Z = 50) and thereby 
should not contribute to an effective deformation. However, the number of neutrons (JV = 58 — 72) is far from closed 
shells (N — 50 or 82), since the shape of the nucleus is due to both kinds of nuclcons wc finally should expect an 
appreciable deformation for these nuclei. Experiment corroborates that because quadrupole-beta values are found to 
be within the interval 0.19 < j3 < 0.29. The most deformed nuclei are the isotopes 118-122 X enon with a quadrupole 
deformation j3 about 0.26 ~ 0.29. These nuclei correspond to the neutron mid-shell region (N = 64 — 66). 
An other common interesting characteristic of these nuclei is the experimental ratio of the energy levels i? C xp = 
E(Ai)/E(2i) which is about 2.3 ~ 2.5 for all these isotopes. This means that these nuclei must belong to the shape- 
phase transitional region between the vibrational limit i? 0X p = 2 and the 7— unstable limit i? cxp = 2.5. 
There are a number of theoretical models that attempt to explain the collective states of the nuclei. The two main 
categories of these approaches are: 

(1) The GBH or geometrical model (see the next section) which treats the even-even nucleus as a quantal liquid 
drop which vibrates and rotates with coupling effects. It is a five dimensional Hamiltonian in which two collective 
variables are devoted to the vibrations of the nucleus and three Euler angles are used to specify its orientation with 
respect to the lab system. Basically it describes the dynamic of the quadrupole deformations of the nucleus but 
sometimes octupole or higher multipole orders are also considered. Some other models such a rovibrational model 
(RVM) can be assimilated to a particular case of the GBH model. 

(2) The interacting boson model (IBM) or algebraic model which assimilates the even-even nucleus to some bosons 
outside the closed shells. These bosons are of type d or s and are made of two nucleons strongly linked by the 
pairing interaction. For odd nuclei the IBM model is replaced by the interacting bosons-fermions approximation. In 
fact there are two versions for this model. In the IB Ml, neutrons and protons are considered as the same bosons 
whereas in the IBM2 these kind of particles are considered as distinct. IBM models can be understood only from the 
group theory and the Lie algebra. Symmetries play a major role in this model. In the IBM1 model the vibrational, 
rotational and 7-soft (or 7— unstable nuclei) limits correspond respectively to the so-called SU(5),U(3) and 5*0(6) 
symmetries. Unlike the GBH, the IBM has a good number of particles (bosons). 

All these approaches contain a certain number of parameters which can be considered as free and therefore can 
be fitted to experimental data to obtain the best possible results. But these parameters can also be derived from 
a microscopic theory. This is a reason why it has no sense to compare models in which parameters are fitted to 
experiment data with others that are based on pure microscopic approaches. 

In our case the GBH model contains seven parameters or more exactly seven functions. Apart from the macroscopic- 
microscopic method which is used to obtain the collective potential energy, the other six functions are evaluated on 
the basis of p ure microscopic models. Similar calculations (but not exactly the same) have already been made in the 
past [l5|-[l6j]. 



II. THE GBH MODEL 



Historically the Bohr Hamiltonian was established as a phenomenological model to interpret harmonic vibrational 
and rotational spectra of the nuclei. Nowadays, the term "Bohr Hamiltonian" is commonly attributed to several 
similar collective hamiltonians. In this respect, we cite Ref. [lij]: "The present-day notion of the Bohr Hamiltonian is 
not very precise. It encompasses a large class of Hamiltonians of which the original Bohr Hamiltonian is only a very 
special case [l5j]-[28[. Here, the GBH means a generic second-order differential Hermitian operator in the Hilbert space 
of functions of quadrupole coordinates. It is the most general collective Hamiltonian using the quadrupole coordinates. 



Making some natural assumptions [14(, it is able to treat large amplitude collective motion. The main advantage of 



the GBH over the IBM model is coming from the fact that it can be derived from a microscopic theory. Two methods 
are usually used to this end: (1) The adiabatic time dependent Hartree-Fock-Bogoliubov method (ATDHFB) which 
leads to quantize with some ambiguities a classical Hamiltonian and (2) The Gaussian overlap approximation method 
associated with the generator coordinates method (GOA+GCM) which gives straightforwardly the quantum collective 
Hamiltonian. Both methods can be applied to different microscopic models such as the mean fields models based on 
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the Nilsson, Woods-Saxon potentials or even self-consistent calculations with Skyrme or Gogny effective interaction. 
In this work we deal with the Generalized Bohr Hamiltonian defined as a sum of three operators: 



H co} = T mb (f3, 7) + T rot (J3, 7, 0i, 9 2 , 3 ) + tfcoiO?, 7) 
where the kinetic vibrational energy and the kinetic rotational energy are given by [l4j: 



T " &(/3 ' 7) = -2^b { ^ 
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T rot (/3,7,0) = 2Lfci.2.3 ^ 

Whereas the collective potential energy ?7 co i of the nucleus is defined as the potential energy of deformation of the 
nucleus (see the following section): 

r and w are given by: w — BppB 11 — Hi , r = SJi^Qfe 

As already mentioned the GBH is a five dimentional Hamiltonian and hence contains five collective variables 
(/?, 7, 0i, 02, 03). The GBH includes seven functions: The collective energy of deformation U co \(f3,j), the three mass 
parameters Bpp, Bp*, _B~ 7 and the three moments of inertia 3i, 92,^3 with respect to the principal axes. All these 
functions are deformation dependent. 
The eigenvalues problem of the GBH has usual the form: 

tfcol^col = ^col^col (1) 

Analytical solutions of the Bohr Hamiltonian can be found in some remarkable cases of potentials: (1) the gamma 
unstable nuclei, (2) the harmonic oscillator potential, (3) the symmetric rotor model. These cases correspond respec- 
tively to the three symmetries SO(6), U(5), and SU(3) of the IBM approach. Obviously these extremes situations are 
ideal cases not met in realistic situations. Shape-phase transitional nuclei occur then between regions of these three 
limiting cases. It is well know that the GBH works better for regions that are far away from closed shells. Because 
in the general case there is no analytical solution, one needs then to solve numerically the Bohr Hamiltonian. Among 
the numerous methods we have chose the one of Libert (with its FORTRAN code) Ref . [24j . The calculations are 
done in two steps: (1) One builds the representative matrix of the Bohr Hamiltonian with the help of a suitable basis 
and then (2) One diagonalizes this Matrix. Beside the collective states obtained with this method it is also possible 
to deduce other observables such as electric or magnetic transitions probabilities, equilibrium shapes for ground state, 
etc... 

III. THE MICROSCOPIC MODEL, SOME NUMERICAL DETAILS: 

In this work, single-particles energies and wave functions are obtained by the diagonalization of the Schrodinger 
equation of the stationary states. 

h4>i(r) = (r) (2) 

€i and <f>i (r) , are respectively, the eigenenergies and the eigenfunctions of the single-particle Hamiltonian h 

h=- (h 2 /2m) VM e// (r)V + V(r) - (n/H) (VW so (r) x p) a + e$ Cotli (r) (3) 

This Hamiltonian contains four contributions. They are respectively: (1) the Kinetic energy operator, (2) the deformed 
central mean field, (3) the spin-orbit contribution, and (4) the Coulomb energy for the protons. The quantities 
M e ff(r), V(r), W so (r), $ CouZ (r) are respectively the effective mass field, the deformed central mean field, the deformed 
spin-orbit mean field and the Coulomb mean field (a denotes here the Pauli spin-matrices). For simplicity we have 
took M(r) = I. 

Single-particle hamiltonian given by ([3]) has exactly the same structure as that obtained from the Hartree-Fock 
self-consistent method with an effective nucleon-nucleon interaction of the Skyrme III type. In our method we have 
simply replaced the self-consistent one body potentials V(r), W(r) by phenomenological deformed mean fields of the 
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neutrons protons 



V = -45.99MeV 
a v = 0.70 fm 
Rv = 6.61 fm 


potential depth 
potential diffuseness 
potential radius 


Vo = -53.22 MeV 
a v = 0.70 fm 
Rv = 6.25 fm 


k = 17.74 MeV fm 2 
a so = 0.70 fm 
Rso = 6.42/m 


spin-orbit coupling 
spin-orbit diffuseness 
spin-orbit radius 


k = 21.13 MeV fm 2 
a so = 0.70 fm 
R so = 5.88 fm 




charge radius 


R ch = 6.25 fm 



TABLE I: Parameters of the Woods-Saxon potential. 



Woods-Saxon type. However for the protons the Coulomb potential has been approximated by the one of a continuous 
liquide drop model with a sharp nucleus surface. For the Woods-Saxon potential the universal parameters' of Ref. 



27| have been used (see also Ref. |28[). This set of parameters claims to be able to reproduce the correct sequences 
of the single-particle levels and also the nuclear equilibrium deformations throughout the entire chart of nuclei [[23 , 
30]]. This set is the one of g| 8 Xe64, it is given in Table|U 

Single-particle states in Eq. @ are solved by the FORTRAN code of Ref. [25j . With these single- particle states we 
calculate the potential energy surfaces of the nucleus (deformation energy) by means of the macroscopic-microscopic 
method. The shell correction is calculated by a semiclassical approach of the Strutinsky method [31(. This means 
that we simply replace the Strutinsky level density by the semiclassical one. This allows us to avoid the well- 
known drawbacks (smoothing parameters) of this method. The pairing interaction is taken into account by the BCS 
approximation. 

We give below some "technical" details of the calculations. 

The deformation energy is defined as the liquid drop energy plus shell and pairing corrections: 

E def (N, Z, 0, 7 ) = E LD (N, Z, p, 7 ) + 5E SC (N, p, 7 ) + 8E SC (Z, p, 7 ) + 5P(N, P, 7 ) + 5P(Z, p, 7 ) (4) 
Where the liquid drop energy is given by: 



3 p 2 Z 2 
E LD (N,Z,p, 7 ) = ~ 



5 roA 1 / 3 



-Ac[^(^, 7 )-l] + [Bc(^7)-l] 



(5) 



in which where B s , and B c are the (normalized) surface and Coulomb contributions to the liquid drop. In the liquid 
drop model we have taken as in Ref. r « 1.275 fm, and C = 52.8(1 - 2.84/ 2 ), I = (N - Z)/(N + Z) 
The microscopic shell corrections SE SC are evaluated separately for neutrons and protons. They are defined as the 
difference between a sharp sum and a smoothed sum (between brackets) of the energy levels: 

SE shell (N or Z, (3, 7 ) = 2 e k - 2 (j^ ) (6) 

with = J £ 9sc{e)de (7) 

— oo 

Here g sc is the semiclassical level density (which by definition does not contain shell effects ) and A sc is the corre- 
sponding Fermi level which is fixed by a constraint on the particle-number. The factor 2 is due to the time-reversal 
symmetry. At last, it is worth to mention that the deformation dependence of the shell correction is contained through 
the eigenvalues tk which depend themselves on the nuclear quadrupole deformation (/3, 7 ). 

The pairing correction SP is evaluated (also separately for protons and neutrons) with the same method as in Ref. [33] . 

SPpairing (N OT Z, p, 7) = P - P (8) 

where the pairing energy P and the " average or smooth" pairing energy P are given by: 



N P 



P = Y, 2 vle k ~ 



k=l 



A 2 

G 



N p /2 

fe=i 



P = -\g sc {\W 



(9) 
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The number of levels used in calculations is defined in formula V3 of Ref . [33| . 

The microscopic moments of inertia and mass parameters are calculated in the usual cranking approximation of 
Inglis-Belyaev: 



^ k = 2h^M^)\ 2 iUvV l + U ^ ] \ k = 1,2,3 (10) 
Da = 2h}^ H — — H -f— TT3 . «,J = /?or 7 (11) 

Here i, j represent the quadrupole deformation parameters f3 or 7, jfe is the single-particle angular momentum and h the 
single-particle Hamiltonian,. The quantity E v = y (e„ — A) 2 + A 2 represents as usual the energy of the quasiparticle 
and u 2 = 1 — v"l is related to the occupation probability of the level v. 



IV. RESULTS: 



From numerous calculations based on microscopic evaluations of the inertial functions (i.e. with no free parameters 
in the GBH), it turns out that generally real difficulties are encoutered to reproduce properly the experimental 
collective levels. Among these difficulties is the fact that the scales of theoretical collective spectra are generally 
too stretched compared to the experimental ones. "Fine structure" such as the order of the levels and their relative 
positions become then a challenge and it is not easy to correct these defects. This is because the seven functions of the 
GBH cannot be deduced directly from experiment and hence cannot be known without ambiguity. Moreover it often 
happens for the same model that good predictions in theoretical spectra are not "corroborated" by good probability 
transitions and vice- versa. 

In the following our theoretical spectra will be compared to the ones deduced experimentally and compiled in the 
websites: http : / /www. nndc.bnl.gov/endsf. 
Our main results can be summarized as follows: 

The spectra obtained by our calculations are given in figures [T] and [5] at the left hand side for each isope. The 
theoretical low energy levels are grouped in sets. 
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FIG. 1: Theoretical collective levels calculated by the General Bohr Hamiltonian with potential energy surface and inertials 
functions evaluated by means of microscopic method using the Woods-Saxon potential. We considere here only quadrupole 
collective level. The parity of the theoretical levels is therefore positive. This is the reason why it is simply omitted. 
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FIG. 2: Continution of Fig. 1 



They are characterized by a structure which is very close to that of the five dimensional (anharmonic) vibrator . 
The triplets (0+,2+,4+) and the quintuplets (0+, 2+, 3+, 4+, 6+) are as a rule, present for all the isotopes. 112_126 le. 
However, experimental data given at the right hand side for each isotope show that only two nuclei belong to this 
type, namely lls Xe and 120 Xe which are well described here . All other nuclei are characterized by doublets (4 + , 2 + ). 
This suggests that these nuclei have a spectrum structure close to the one of the 7 unstable nuclei of Wilets-Jean 
model which predicts the doublet (4 + , 2 + ). 

We recall that the Wilets-Jean model is essentially based on the hypothesis that the potential energy does not depend 
on the axial assymmetry 7 [35[], i.e. dU co i(/3 , / dj = with a strong minimum out of the spherical shape (i.e. the 
minimum occurs for (3 ^ 0). These nuclei are oten called as to be 7— soft. In this respect, it is possible to consider 
the harmonic vibrator as a particular case of the Wilets-Jean model. In effect, the Wilets-Jean model is defined by 
the condition dU co i((3 , 7) / — which is filled by the harmonic vibrator. The only difference comes from the fact 
that the minimum of the collective potential energy occurs for the spherical shape whereas in the W-J model it lies 
elsewhere. Consequently, it is then not surprising that among them two are of vibrational type. 
It is easy to explain why our theoretical values are close to the ones of the anharmonic vibrator: 

(i) The potential enegy depends very little on the axial assymmetry parameter 7 (as for Wilets-Jean model) 

(ii) The deformation energy (defined as the difference between the energy for the spherical shape and the one obtained 
for the equilibrium deformation) is too weak. In other words these nuclei are theoretically too "soft" in the (3 degree 
of freedom and the well is insufficiently pronounced for f3 ^ to obtain Willets-Jean potential type. 

(ii) Moreover the mass parameters Bpp, Bp^, B 77 vary very little in the vicinity of the minimum. Therefore as in the 
original Bohr model the mass parameters can be considered as constant. 

Above the doublets we find either a set of the type (2 + , 3 + , 4 + , 6 + ) or of the type (0 + , 2 + , 3 + , 4 + , 6 + ). The experimental 
energy ratios E(4^)/E(2^) lie between the the 7— soft and the vibrational limits in Fig|3]but are much more closer 
to the 7— soft limit especially for the last isotopes. The theoretical values are not so far away from the experimental 
ones but their behaviours as a function of A seem to be opposite from each other. 

The first satisfaction of these calculations is the correct scale of the spectra without any kind of correction. In other 
words the usual defect of the stretching [l4| of the spectra is not present. This seems essentially due to the Woods- 
Saxon mean field which gives the right values for the mass parameters. In earlier similar works, for example the one 
of Ref.[l5| the stretching of the collective spectra is "cured" by introducing the pairing vibrations in the calculations 
whereas in Ref. [l6[ it is corrected by reducing artificially the pairing strenghs by 20%. The second satisfaction is 
the calculated levels 2^, and the doublets (4+, 2^) which are for the most nuclei very close to the experimental ones. 
Moreover as already noted the energy ratios given by our model in Fig|3]are quite close to the experimental ones. 
Nevertheless we obtain also bad results such as for the mean values of the ground states quadrupole deformation 
(/3) in Fig 01 The experimental values are deduced from transitions probabilities B(E2; + — > 2 + ) from the website 
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FIG. 3: Energy ratio )/E(2+) for the family of isotopes \f- 126 Xe. 



http : I /www. nndc.bnl.gov/be2 whereas the deformations given by the FORTRAN code come from the collective 
wave function (B(E2; + — > 2+) values are not calculated by the present version of the code). 




FIG. 4: Root mean Square value y/< /3 2 >of the quadrupole deformation in the ground state as function of the mass number 
for the Xenon isotopes. 



CONCLUSION 



Microscopic calculations based on the Woods-Saxon mean field were performed to find the potential energy and 
the six inertial functions entering into the collective Bohr hamiltonian. We have considered the light and medium 
Xenon isotopes of the region 112 < A < 126. Then, the Bohr hamiltonian has been diagonalizcd without any fitting 
of parameters. The resulting spectra were found to be quite close to the experimental ones. Moreover the values 
of the mass parameters seemed to be correct in magnitude giving a good scaling in the collective spectra for this 
region. Most of theoretical collective levels E(2f) were found to be in good agreement with the experimental ones. 
The theoretical energy ratios E(4^)/E(2f) were also fairly well. However, the major difference comes from the fact 
that our spectra belong to the (anharmonic) vibrator type whereas most of experimental spectra exhibit doublets of 
the Wilets-Jean model. Both alike are of the same type (7— soft), however contrarily to the Wilets-Jean model, our 
collective potential energy does not posses a sharp minimum for a deformed shape (/? ^ 0). In other words the shell 
correction is not sufficiently strong to modify significantly the potential energy of the liquid drop model. Consequently 
some further studies seem to be necessary in order to correct this defect. 
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